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' We study the two-dimensional Ginzburg-Landau model of a neutral superfluid 

. in the vicinity of the vortex unbinding transition. The model is mapped onto an 

| effective interacting vortex gas by a systematic perturbative elimination of all fluctu- 

ating degrees of freedom (amplitude and phase of the order parameter field) except 
the vortex positions. In the Coulomb gas descriptions derived previously in the lit- 
erature, thermal amplitude fluctuations were neglected altogether. We argue that, 
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q | if one includes the latter, the vortices still form a two- dimensional Coulomb gas, 



but the vortex fugacity can be substantially raised. Under the assumption that 
Minnhagen's generic phase diagram of the two- dimensional Coulomb gas is correct, 
our results then point to a first order transition rather than a Kosterlitz-Thouless 
transition, provided the Ginzburg-Landau correlation length is large enough in units 
of a microscopic cutoff length for fluctuations. The experimental relevance of these 
results is briefly discussed. 
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1 Introduction 



The critical behavior of two-dimensional (2D) systems with a continuous internal sym- 
metry has been a most puzzling problem for a long time. Simple physical realizations are 
superfluid or superconducting thin films, which on a phenomenological level can both be 
described by a complex order parameter field the boson or Cooper pair "condensate 
wavefunction" . The internal symmetry is then a U(l) gauge symmetry acting on the 
phase of 

Early theoretical work showed that the "usual" criterion for superfluid order in 3D 
bulk systems JTJ, namely long-range order (LRO) of the field is not satisfied in 2D 
films at any finite temperature @. The reason are low-energy phase fluctuations, the 
Goldstone modes related to the broken gauge symmetry, leading to \& correlations which 
decay algebraically to zero with distance. However, it was quickly realized || that this 
"quasi" -LRO was sufficient to insure superfluidity. The correct criterion for superfluidity 
turned out to be rather a nonvanishing stiffness with respect to long-wavelength phase 
fluctuations (the helicity modulus T 0) than true LRO in 

As a matter of fact, experimentally even very thin 4 He films of a fraction of an atomic 
layer showed clear signatures of superfluidity (see refs. in ||). 2D superconductivity was 
predicted |], [7|] for sufficiently "dirty" samples (high normal sheet resistance) and observed 
in continuous and granular thin films (see e.g. || |{| [H| [TT]|). 

The nature of the transition from this superfluid phase to a high- temperature phase 
with exponentially decaying \1/ correlations was clarified by Berezinskii W2\ and Koster- 
litz and Thouless (KT) |13|, [14]]. They realized the essential role of vortices, i.e. phase 
singularities with nonvanishing winding number ("vorticity"). The vortices interact log- 
arithmically at large distances, thus at low temperatures they appear as bound pairs of 
zero total vorticity only which do not change the algebraic decay of the \1/ correlations. 
However, at some finite temperature T v , the largest vortex pairs in the system start to 
dissociate (unbind) by a collective screening mechanism, and the free vortices lead to an 
exponential decay of the correlations. (We use the notation T v here in order to avoid 
confusion with a "conventional" critical point which will appear later.) 

Usual phenomenological models of superconducting or superfluid films are the Ginzburg- 
Landau (GL) model or the XY ("planar rotator") model |7|, [15[ , respectively (see also 
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section 2). The GL model is generally believed to provide an appropriate description of 
superconductors close to the bulk transition since there is a microscopic derivation in the 
framework of Gor'kov theory (see e.g. refs. [|HJ). The XY model may be considered as the 
"phase only" limit of the GL model, in which the modulus of \I/ is fixed and only its phase 
is allowed to fluctuate fll7| . Together with the idea of thermally excited vortex loops, it 
has been succesfully employed to describe the superfluid transition of 4 He in 3D [TS) and 
2D 0, |15ll ) but there is no comparable microscopic justification as for the GL model in 
the case of superconductors. 

Since the KT transition is essentially caused by the interacting vortex system only, it 
appears in its purest form in a corresponding model system of pointlike particles in 2D 
with a logarithmic interaction, the 2D Coulomb gas (2DCG) model (for a review see e.g. 
||T9| ). This model depends on two independent, dimensionless parameters, a temperature 
T CG and a fugacity z CG of the particles. The GL and XY models may then be thought of 



as particular realizations of the CG model, represented by particular z CG (T CG ) lines [19]. 

Based on the ideas of Kosterlitz and Thouless, a systematic renormalization group 
(RG) theory of the KT transition has been developed, usually starting from an XY or CG 
type description |2(], |21fl . In the CG picture, the KT RG equations are flow equations in 
the z CG - T CG plane; the z CG (T CG ) line of a particular realization serves as initial condition 



of the RG flow. In the RG framework it could for instance be shown |22fl that the helicity 
modulus has a finite value T(T — * Tj~>) ^ just below the transition and then drops 
discontinuously to zero. Furthermore, the ratio 27rT(T)//cB^ tends to the universal value 
4 at the transition (T — > T^~>). This result may also be expressed in terms of pure CG 
quantities: general hydrodynamic arguments [23| show that the quantity (eo(T CG )T CG ) _1 
(e being the k — > limit of the CG dielectric function) corresponds exactly to the above- 
mentioned ratio 27rT (T) / k^T at all temperatures, in particular its discontinuity is the 
same. 

This famous "universal jump" prediction is one of the central results of the KT RG 
theory and has been verified to an impressive degree by measurements on 4 He films || 



as well as by numerical work on XY type models (see e.g. [23] and refs. therein, refs. 



in ||25|| ) and CG type models ||26|| . Quite convincing evidence for KT universal behavior 
has also been obtained for artificial superconductor networks, more specifically weakly 
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coupled Josephson junction arrays (JJA's) |27[] and wire networks [25], which are both 
well described by XY type models (see section 2). 

The situation for continuous superconducting films appears to be less clear. Experi- 
ments on high-T c films (YBCO fLO), (TTJ] ) and dirty conventional films (Al [§, In/InO H) 
films provide rather good evidence for a discontinuous jump of the helicity modulus T at 
the transition, however in both cases they seem to point to slightly larger values of this 
jump. (We note that in one case — ref. |T(| — the authors explain the observed behavior 
of T in a completely different way, namely by percolation effects due to the granularity 
of the film.) 

Unfortunately, in the literature apparently no numerical results are available on the 
critical behavior of the 2D GL model. Usually, the 2D GL model is assumed |T7| , |T5H 
to have the same universal properties as the XY and CG models, since small amplitude 



fluctuations of the order parameter \l/ have been shown |29| to be irrelevant in the RG 
sense. These fluctuations are then neglected altogether, and one is left with a gas of GL 



vortices (termed "Ginzburg-Landau Coulomb gas" (GLCG) in the literature pUL [19]]; we 



will call it the "bare" GLCG since fluctuations are not yet included) whose interaction 
is logarithmic at large distances and cut off roughly at the GL correlation length £. 
However, it is easily seen (see subsection 2.1) that right at the expected vortex unbinding 
transition, the amplitude fluctuations cannot be considered as weak anymore and thus 
may well change the critical behavior. 

A further, conceptual problem of the "bare" GLCG is that the phase space division A 
of a vortex is not well defined without additional arguments ||19|| . To get a rough idea of 
the size of A, one may argue that neglecting fluctuations at length scales < £ corresponds 
to a lattice regularization with lattice spacing ~ £. Consequently, the vortex phase space 
division is expected to be ~ £ 2 , since a vortex can be located at any of the plaquettes of 
the lattice. A more systematic approach to this question will be discussed at the end of 
this paper, in subsection 4.4. 

The above qualitative argument implicitly assumes that the only length scale which 
enters the problem is £, in particular the z CG (T CG ) line should be independent of any 
microscopic cutoff length for fluctuations. In contrast to this, in the present paper we 
study precisely the influence of small- wavelength (<£) fluctuations on the transition. 
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In doing so, we still assume that a description of the transition in terms of an effective 
interacting vortex gas is justified. We then take all thermal fluctuations systematically 
into account by a perturbative elimination of all degrees of freedom except the vortex 
positions. Our results will indicate that these fluctuations strongly increase the density 
of vortices at the transition or, in other words, they shift up the z CG (T CG ) line. 

Or, already KT in their original paper [H noted that their approximations are justified 
only when the vortex system is dilute enough at the transition, i.e. when z GG = z CG (T GG ) 
is small enough. Outside the region of small z GG , universality will hold only as far as 
the RG flow remains KT like in a topological sense. On the other hand, Minnhagen 
|3T| (see also [^9, |32], |33|) investigated the whole phase diagram of the 2D CG model by 
an extension to larger z CG of the KT- Young self-consistent screening procedure |14], |34j] 
of deriving the RG equations. Minnhagen's modified RG equations are nonholonomous 
integro-differential equations which in the small z CG limit reduce to the KT equations. 
Moreover, for small enough z CG his corrections to the latter are irrelevant in the RG 
sense, i.e. Minnhagen's equations reproduce the same universal properties in this region. 
At larger values of z 0G however, his equations predict a first order transition line which 
ends at some critical point (z CG ,T GG ) w (0.029,0.204). The KT line joins the first order 
line smoothly from below at some point (z CG ,T GG ) w (0.054,0.144) (fig. [I], qualitatively 



adapted from f33j). The superfluid transition is KT like up to this point and first order 
further above. 

The first order part probably may still be interpreted as a vortex unbinding transition, 
which in contrast to the KT transition involves the simultaneous dissociation of a finite 
fraction of all bound vortex pairs in the system. Its characteristics are nonuniversal, for 
instance the jump of the above- mentioned quantity (eo(T CG )T CG ) _1 depends on T GG in 
the first order part of the transition, and it is larger than the KT value 4. 

In fig. |l|, we have included the z CG (T CG ) lines representing the XY model and the 
"bare" GLCG, where for the latter we used the functional form proposed in (see also 
subsection 2.3). The XY line definitely lies in the KT regime, there is thus no contradiction 
between the Minnhagen theory and the numerically well established fact (see references in 



24| and P5| ) that the XY model displays the universal properties of a KT transition. The 
GLCG line intersects both the KT and the first order lines, but the superfluid transition 
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is still KT like. Our aim here is to argue that inclusion of fluctuations may shift the 
z CG (T CG ) line even further up in the first order regime. 

However, the Minnhagen scenario must still be considered as somewhat speculative. 



Some effort has been made in the past (see |25|, ^] and references therein) to decide 
whether a first order transition exists in a modified XY model with a "truncated" cosine 
interaction between neighboring phase angles (which presumably also corresponds to a 
higher z CG (T CG ) line than the XY one in fig. [I]) but the conclusions obtained by different 
authors are contradictory. Perhaps the strongest case in favour of the Minnhagen scenario 
could be made by Jonsson et al. P5| . By a finite-size scaling analysis, they established 
the existence of a critical point above the KT line, which they identify with (z GG ,T GG ). 

In summary, in this paper we try to treat the mapping of the 2D GL model on the 
corresponding interacting vortex gas in a somewhat more complete and systematic fashion 
than it has been done so far in the literature. Thereby we find that fluctuations strongly 
enhance the vortex fugacity at the transition, which may drive the transition first order, 
provided the Minnhagen szenario is correct. 

The paper is organized as follows: in section 2, we start with some general remarks 
about the GL description of superconducting films and networks, and also define our 
modified notion of the "Ginzburg-Landau vortex gas" (GLVG). The details of the main 
calculation are contained in sections 3 and 4. The idea is to first investigate some kind 
of "saddle point" configuration of the field for given vortex positions (section 3), which 
corresponds roughly to the "bare" GLCG discussed previously in the literature. In a 
second step, we include thermal fluctuations around this configuration in a gaussian ap- 
proximation (section 4) and derive, as our central result, the z CG (T CG ) relation for the 
GLVG. Section 5 contains a short summary and conclusions. 



2 The effective Ginzburg-Landau vortex gas 



6 



2.1 Ginzburg-Landau description of superconducting films 

The starting point of this paper is the GL free energy functional of a complex order 
parameter field ^ in 2D, 



H[V\ = |d 2 r ja|^| 2 + ^| 4 + 7 |W| 2 j . (1) 

A physical system can be considered as 2D, if the thickness of the film is smaller than 
some minimum wave length of ^ fluctuations, such that variations of \I/ perpendicularly 
to the plane are negligible. Note that we also omitted any coupling to the magnetic vector 
potential A, i.e. ([I]) actually describes a neutral superfluid. Nevertheless, as argued in 
the introduction such a local free energy is probably more appropriate for superconduct- 
ing films than for 4 He films. At the end of this section, we will indicate under which 
circumstances the omission of a coupling to A can be justified. 

7i is used to define statistical-mechanical quantities like the partition sum 

Z oc J D* e -W*BT (2) 

The coefficients a, /3, 7 in TC are in general temperature dependent quantities: as usual, we 
assume that a vanishes at some "mean field" critical temperature T c0 , and that a < for 
T < T c o- Furthermore, for (Q) to be defined, necessarily /3, 7 > 0. There is some freedom 
in the interpretation of the parameters a, f3, 7 and the order parameter field usually, 
one chooses 7 = h 2 /2mn where rnn is the effective in-plane mass of the carriers (electrons 
or holes) and interpretes |\l/| 2 as the local density of carriers in the condensate. If we 
disregard boundaries, (|I|) assumes its minimum for the homogeneous field configuration 



= y\a\/P and for any other configuration which differs from by an arbitrary, 
spatially constant phase factor. 

For later convenience and clarity we introduce a dimensionless order parameter field 



ip := J/3/\a\ty normalized such that ^ = 1. With 



and 



f==JS (3) 



1 /_ , ,o\2 



m= I d 2 r j— (1 - H 2 ) +|V^| 2 | (4) 
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we can then rewrite the exponent of @ as 

= jHty] + const. . (5) 

"const." stands for a ^-independent term which is irrelevant for the thermodynamics. £ 
is the GL "correlation" length which measures the length scale on which ip relaxes to 
ipoo = 1 away from a perturbation. We will see later that it is also the correlation length 
for thermal fluctuations of the amplitude K measures the stiffness of the ip field over 
temperature and its inverse will later play the role of an effective statistical-mechanical 
temperature. Note that 1/K and £ both diverge as T approaches T c0 from below, since 
\a\ — > 0. Both effects tend to enhance fluctuations, so if the KT transition is not preceded 
by some other transition to a disordered high-temperature phase, there must be a vortex 
unbinding transition at some temperature T v below T c0 . 

Expression (|J) looks as if £ were the only length scale in the problem and could be 
eliminated by a simple rescaling of all lengths, i.e. of r. This is in fact the case in usual 



GL theory [16|. However, in order to obtain well-defined, finite results for statistical- 
mechanical quantities like the partition sum (Q) or ip correlation functions, we have to 
limit the number of degrees of freedom by a UV regularization which introduces a second 
length scale. For definiteness, we may for instance put the model on a square lattice 
with spacing a, 

m = i- 2 E (i - i^i 2 ) 2 + E - n 2 , (6) 

K i (ij) 

where denotes a sum over all pairs of nearest neighbor sites. Alternatively, we may 
supplement the continuum version ([|) by the corresponding momentum cutoff prescription 
(i.e., all momenta in the Fourier representation of (^) restricted to the first Brillouin zone 
[— -, ^] 2 ). We will assume that both prescriptions are essentially equivalent and use them 
in parallel. 

Expression (|6|) clearly shows that £/a is a second (besides K) dimensionless parameter 
which enters the model. In the limit ( < a, the first term in ([]) suppresses amplitude 
fluctuations of the order parameter away from the value \ipi\ = 1 and (H) reduces to the 
XY hamiltonian, 

#xyM = E \^ - ^ I 2 = 2 E (i - cos (^ - . (7) 



s 



where 0j is the phase of ipi. In this respect, the XY model is the "phase only" £ — > 
limit of the GL model. As mentioned in the introduction, its predictions agree well with 
measurements on Josephson junction arrays (JJA's) and superconducting wire networks. 
From our point of view, the reason for this success is that in these artificial networks the 
spacing of the underlying lattices provide a macroscopic cutoff length a which can be tuned 
independently of £. Weakly coupled networks P7|, which are preferred experimentally 
since their transition temperature is well seperated from the bulk transition then always 
lie in the regime £ <C a. Note that in a model of type (||) for weakly coupled JJA's, 
£ is not equal to the bulk GL correlation length but usually much smaller (as a matter 
of fact, £ 2 /a 2 is essentially proportional to the Josephson coupling energy between two 
grains over the condensate energy of a grain). Consequently, there is no contradiction to 
the fact that usually a is smaller than the bulk £ in these systems. 
The main points we want to stress now are the following: 

1. In continuous superconducting films (as opposed to networks), £ is not independent 
of a. In fact, it always obeys the inequality £ > a. 

2. If £ > a, thermal amplitude fluctuations are quite strong close to the transition. 

In the remainder of this paper, we will then analyse in detail the effect of these fluctuations 
and discuss in which way they may change the critical behavior. 

In order to demonstrate point 2., we disregard for the moment the phase degrees of 
freedom and estimate the local fluctuations of the amplitude in a gaussian approxima- 
tion. An expansion of the potential term in the hamiltonian (^) shows that amplitude 
fluctuations 5\ip\ = 1^1 ~~ 1 have a "mass" 2/£ 2 . Therefore their mean square value at 
some "temperature" 1/K may be estimated as 

((5|y>|) 2 ) _ j_ /■*/« d 2 k i 

(H) 2 ~ 2K J-n/a (2tt) 2 2/e + k 2 
1 r^/* 2 d{k 2 ) 1 
~ 8^K Jo 2/£ 2 + k 2 ~ 

where in the second line we just have approximated the quadratic Brilluoin zone by a 
circular one of the same area (27r/a) 2 . At the presumed vortex unbinding transition, 
one expects the dimensionless "inverse temperature" K to be roughly of the order of 1, 
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which implies that for £ > a the fluctuations are of the same order of magnitude as its 
expectation value squared, {\ip\) 2 ~ 1. 

To show that for superconducting films always £ > a (point 1. above) we first have to 
understand the meaning of cutoff a in this context. In fact, the Gor'kov derivation (f37|, 



see also the reviews |16[]) of the GL functional (JT|) from BCS theory does not immediately 
yield a local form like ([I]), but involves integral kernels whose range plays the role of the 
cutoff a. Both a and the GL correlation length £ can be calculated in this framework and 
one obtains roughly 

i« , 1 . (9) 

where \ pa (1 + h^nk-sTr) -1 is a number < 1 for dirty superconductors and = 1 in the 
clean limit |16|| . An immediate consequence of (^) is the validity of the asserted relation 
£ > a at any temperature T. £/a diverges as T approaches T c o from below; note that 
one usually even assumes £ ^> a to justify a local GL description. From BCS theory we 
can also estimate the value of £/a in the interesting region close to the presumed vortex 
unbinding transition: using fl3|), (H) one can show that 



h 2 nor, 

210, 



n 2 B [A 2 ] m c 



(10) 



a ~ ^ 2m||A;BTfs: ~ ^ T[°K] X my ' 

where n2o[A 2 ] is the number of carriers in the film per A 2 and mil their in-plane effective 
mass. At the transition (T = T v ), we assume again K pa 1 and fliCf ) yields a £(T V ) 
which is appreciably larger than a unless the film has at the same time a high T c o, low 
carrier density ri2D and high in- plane carrier mass m\\ (note that these conditions may 
be realizable in films of high-T c material whose thickness is a few unit cells). 

The last point we want to briefly address in this subsection is the omission of a 
vector potential term in ([I]). A priori, such a coupling is important also in the absence 
of an external magnetic field, since thermal ip fluctuations are accompanied by local 
supercurrents and therefore generate local magnetic fields. 

The argument for the conventional scenario (no fluctuations except vortices) is well 
known || |7| and goes as follows. It was shown by Pearl |38[| that the logarithmic vortex- 
vortex interaction in superconducting films is magnetically screened at a length 
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where Al is the bulk London penetration depth and d the film thickness. For this reason, 
KT originally argued that a vortex unbinding transition should not occur in supercon- 
ducting films since the logarithmic interaction at arbitrarily large distances is essential for 
the mechanism. However, it was noticed later || 0] that for small d and sufficiently close 
to T c0 (note that Al diverges as T — > T c0 ), A s can become of the order of millimeters, i.e. 
comparable to typical sample sizes. Screening then should not wash out the transition 
to a larger extent than effects of finite system size and is therefore disregarded. The 
above authors restricted their argument to dirty films with a large Al(T = 0) but one can 
easily convince oneself using the GL relationship A£ = h 2 c 2 p/32n\ a|7e 2 that for any film, 
independently of its specifications, 

(hc/e) 2 1.2 cm 
s = 8irKk B T ~ T[°K] K ■ [ ' 

At T = T v , A s should thus always be roughly of the order of 1 mm. 

In this paper, we consider — in addition to vortices — the effect of amplitude and 
phase modulations which are of course also accompanied by supercurrents. Screening 
effects lead in this case to an additional "mass" term ~ 2/A;? for fluctuations, in particular 
the massless phase (Goldstone) modes of the neutral superfluid become slightly massive. 
However, since our calculations will involve fluctuations of wavelength < £ only and since 
we can assume £ A s because of (0), we believe that we may disregard this effect, too. 
A more detailed discussion of this assumption may be interesting, but we will not further 
enter this question here. 

2.2 Elimination of short-wavelength fluctuations: mapping on 
the two-dimensional Coulomb gas 

In this section, we outline the basic ideas of our approach. In particular, we try to clarify 
how the mapping of the 2DGL model at finite T on its associated vortex gas can be 
carried out in a systematic way. 

Our aim is to study the presumed vortex unbinding transition which should show up 
in thermodynamic quantities like the partition sum (|2|). To this end we try to transform 
(FJ) into the statistical-mechanical partition sum of an interacting vortex gas, i.e. we try 
to eliminate any degrees of freedom except the vortex positions. 
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On the technical level, the main step is to introduce a new lattice cutoff a > a and 
to define an effective action S[ipo\ on a coarser lattice with constant a (i.e., of the long- 
wavelength field components ipo) by integrating out the short-wavelength components ipi 
which contain wave numbers between n/d and n/a: 

Z<x[ ^Dfoe* 1 , (13) 

a 

e §[i>o] := f 2^ ie -f^+^ . (14) 
J l< k <l 

If a is chosen large enough (a>a,£) then we expect (0) to be dominated by "saddle 
point" configurations ipo of H with vortices at given positions Rj (and with associated 
vorticities mi) as constraints. The effective phase space of the resulting vortex gas will be 
formed by the plaquettes (with area a 2 ) of the coarser lattice which is important for the 
definition of the vortex entropy and thus for the thermodynamics. Non-singular phase 
fluctuations ("spin waves") of wave numbers k < n/a are also still contained in (0), but 
they decouple from the vortices like they do in the pure XY model and play no role in 
the transition, so we will disregard them (see also subsection 4.2). 

On the other hand, if a is not much larger than max(a, £), then the effective action 
S[ipo] may be determined to a reasonable approximation by a gaussian approximation to 
the functional integral in flT4]). a may then be thought of as the smallest length scale on 



which vortices are well distinguished objects. In subsection 4.2 we will realize the infrared 
cutoff 7r /a by a mass 2/£ 2 of the field ip± ("soft" cutoff), where £ is the GL correlation 
length. Both quantities will then turn out to be related by a 2 = a 2 + 27r£ 2 which is in 
good agreement with our general discussion in the introduction and satisfies the above 
criteria. 

Disregarding phase fluctuations of wave numbers k < n/a as stated above, the vortex 



part of (p~3[ ) can be expressed as 



z-s^Efn/^^n, (15) 



N {rrii} \i=l-' fl2 



where the determination of S[i/jq N ^] as a function of the Rj,mj from (|14D is the main 
technical task in this paper. iV is the total number of vortices (of mutual separation > a) 
in the system. We will later see that thermodynamically (i.e., in the partition sum (|l5|)) 
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only vorticities = ±1 are important and, moreover, below T v only "neutral" vortex 
configurations (obeying J2i m i — 0). 

As we will see in subsection 4.3, for vortex distances much larger than a, £[0o ] has 
an asymptotic behavior of the form 



S[4 N) ] ~ ^ ( E ™ t m 3 In |R * a Rjl + iV> CG j , (16) 
i.e. it behaves as a neutral two-dimensional Coulomb gas (2DCG). The latter is generally 



defined by a partition sum of the form [19 



^E^E'fn/^W-L Em^W^-^l + ivH , (17) 

N ly -{m t } \i=l J J 1 \%<j ) 

where J2'{ mi } now i s a restricted sum over all neutral configurations of iV vortices with 
rrii = ±1. A is the "phase space division" of a CG charge, T CG the dimensionless CG 
temperature and — 2/z CG may be interpreted as the creation energy of a neutral pair at 
distance R c . Note that in ([17]) the interaction involves \n + x := max(0, ma;), i.e. it is zero 
at smaller distances than R r . Such a cutoff is essential for (1l7D to be well-defined, but 



it can be realized in different ways ||19|| . It is easy to see that despite the appearence of 
the four parameters T CG , /z CG , A and R c , (|T7J) actually depends only on T CG and on the 
"fugacity" 



z cg = th c M CG /r CG 

A v 



of the CG charges. The values which were used to draw the z CG (T CG ) lines in fig. |I] are 
taken from refs. (XY model) and |TJ5], ("bare" GLCG), and are collected in table 



I. For the effective GLCG defined by ([15]), (|16D we have by definition A = a 2 , but the 
short-distance cutoff is not yet specified in (|16]). We assume here R c ~ a, such that 
simply z CG = e M ° G ^ T ° G . This is reasonable since a is of the order of the vortex core size. 
To demonstrate that the results for z CG (T CG ) are not too sensitive to the choice of R c (nor 
presumably to the precise cutoff procedure), we have changed the values of the cutoff 
distance R c in the case of the XY model and the "bare" GLCG by factors of 2 and 1/2. 
Note that there is an accompanying shift of /i CG : the change to some other cutoff R' c in 
( |TT| ) means that we have to replace \ia + (R/R c ) i— > \n + (R/R' c ) + ln(i?^./_R c ) which (using 
the charge neutrality condition) leads to 

/ „/ \ 2-l/2T CG 

*"^it) ■ (19) 
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The results are shown in fig. || and present no essential changes compared to fig. [!]; the 
effects we will be concerned with later on are much more pronounced. 

3 The "saddle point" configuration 

As we argued in subsection 2.2, the long-wavelength component i/jq of the field is essentially 
of the form of a "saddle point" configuration of H for given vortex positions Rj and 
vorticities m,. So our first step will be a detailed investigation of these configurations and 
in particular of their energies Htyo]. In section 4 we will then proceed calculating the 
effective action S[ipo] of (|T4T), by adding corrections due to thermal fluctuations around 
these saddle points. 

Since in this section we are not considering fluctuations, we can work in the contin- 
uum limit a — > (such that £ is the only length scale of the hamiltonian (01) ) without 
encountering divergencies. Only in subsection 3.5 we will reintroduce a discrete lattice 
and discuss, in which way our results then are modified. 

3.1 Separation of the vortex degrees of freedom 

In order to see how vortices, i.e. singularities of the phase with finite winding numbers 
("vorticities") can be imposed as constraints on the field r/>, we write the latter in terms 
of real fields p,9 ("amplitude" and "phase"), 

$ = pe ie (20) 

and express the hamiltonian (||) in terms of p, 9: 

H\pA = /^|^(l-p 2 ) 2 + |Vp| 2 + p 2 |V^| 2 | . (21) 

The phase gradient can be split into its longitudinal and transverse parts, 

^6 = Vtf - n x V$ , (22) 

where n is the unit vector normal to the plane. is a nonsingular phase field representing 
the "spin waves" , and $ is a "vortex potential" which satisfies the Poisson equation 

V 2 $(r) = -2tt J2 m i S ( r - R i) ( 23 ) 

% 
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with point-like integer "charges" (vorticities) at the positions Rj, corresponding to the 
singularities of the phase field 9. (6 changes by 27rmj upon going counterclockwise around 
Rj). $ is determined by ( p3|) only up to a harmonic function; any harmonic contribution 
can however be absorbed into so that we may choose the particular solution 

$(r) = -^m 4 ln^— 5±1 (24) 

which finally renders the splitting (|2"2] ) unique. In order to make the argument of the 
logarithm in fl24|) dimensionless we have introduced the only length scale £ of the prob- 
lem. Note however that for a "neutral" vortex configuration {J2i m i = 0), $ does not 



depend on £. The physical content of (22]), ([24]) is that we have separated the vortex de- 
grees of freedom (expressed by their positions Rj and vorticities mi) from the remaining 
nonsingular phase configuration 



3.2 The vortex core structure (1-vortex problem) 

Besides the phase field which defines the vortices, our hamiltonian contains the amplitude 
field p which is strongly coupled to 9 close to the vortex centers R^: far away from the 
centers (V8 small) it is expected to have values p « 1, whereas it vanishes right at the 
vortex centers, p(Rj) = 0, because otherwise the "vortex core" energies would diverge 
logarithmically in the continuum limit. We call vortex core the regions of size ~ £ around 
Rj where p significantly differs from 1. Since the structure of the core regions is essential 
for the following, we will investigate here in detail a single isolated vortex of vorticity m 
centered at the origin. This subsection is sort of a summary of results taken from the 
literature which are relevant in our context. 

Because of the isotropy of the problem we employ polar coordinates r, 0; in the saddle- 
point configuration of the hamiltonian (pT|), the phase field is then (up to an additive 
constant which we choose equal to zero) given by 

# m (r,0) = m0 (25) 

and the amplitude p m is a function of r alone. Inserted in (|2l|), this yields a reduced 
hamiltonian of the 1-vortex problem, 

H m [p m ] = J dr2nr |^(1 - p 2 J 2 + + ^-pA . (26) 
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Its minimum solution is determined by a vanishing functional derivative with respect to 

Pmi 

n _ 1 6H m _ d 2 p m 1 dp m m 2 J_ 2 

4vrr5p m (r) ~ dr 2 r dr r 2^+£2^W P»J , l^J 

together with the boundary conditions 

p m (0) = , p m {r) -> 1 for r -> oo . (28) 

The solution p m (r) rises monotonically from at r = to 1 at r — > oo. To determine the 
asymptotic behavior of p m at small and large r, one may proceed as follows: 

1. r -C £: here p m (r) <C 1, so we linearize ( p7|) in p m and insert an ansatz p m (r) ~ 
c m (r/£) a with a > and c m a numerical constant, which yields 

~ ^ + - ^p m ~ c m (a 2 - m 2 )^ . (29) 

Thus a = \m\, and c m is a constant of order 1 which has to be determined numeri- 
cally from the full solution of fl2T|), (|28|). For m — 1, reference |30[ gives 



lim ( ^-^-p 2 m= A = ^ = 0.108 , i.e. c m=1 = 0.582 . (30) 
r->o \ 2irr ar ii 



2. r ^> here w(r) := 1 — p m (r) <C 1, so we linearize ( f27|) in u and insert an ansatz 
u(r) ~ c'^r/O' 3 with (3 < 0: 

d 2 w 1 aw m 2 . . 2 

~ 4(/5 2 -^ 2 )^ + ^-24^. (31) 
For r — > oo, the last two terms are the dominant ones, such that (3 = —2 and 

4 = ™ 2 /2. 

In conclusion, the asymptotic behavior of p m is given by 

c m (r/0 H for r<£ 
p m (r) ~ <( m 2 ' , c m=1 = 0.582 . (32) 

for r»e 



Sketches of p m for m = 1, 2 are given in reference pJO]. For later use we note that the first 



line of fl32|) together with 8 m = mcj) imply that close to the vortex center (r -C £) the field 
0m — Pm£ t6m has the simple power form 

, c m (f/0 |m| for m > 

mK/q) , (33) 

c m (f*/0 for m < 
16 



where f = re 1 ^ is a complex notation for the coordinate in the plane. Note in particular 
that the complex field is perfectly smooth even in the center of a vortex, a singularity 
only appears when one looks at the phase separately. 

Now consider the different contributions to the total energy H m (|26|) . Since for a 
circle-shaped system H m diverges with the radius r c as ~ 2irm 2 lnr c , we split off this 
size-dependent term: 



H m [p m ] ~ 2vr (m 2 E x {m) + E 2 {m) + + m 2 In T -A (34) 



for r c ^> where 



E x {m) := lim I / dr^ - In ^ ] = - f°° dr^ In T - , (35) 



r £ y jo dr £ 

dr 



E 2 (m) := / drr[^] 2 , (36) 
E 3 (m) 



POO 

y o ^-(l-p 2 ,) 2 . (37) 

The second expression for E\{m) follows by partial integration. 

E 3 (m) can be calculated analytically by the following trick (we follow an idea of 
reference PD|): let p m be the solution of (p7j), ( p8|) and let pm )a {r) := Pm{r/a). Then the 
energy H m [pm,a\ is minimum for a = 1 . To yield finite results it must again be regularized 
by a finite system radius r c ^> £. After a variable change r i— > ar it reads: 

H m [ Pm , a }= r rc dr27rr\^(l-pir+(^)\^pl) . (3* 



2 ^2 v rmJ ' \ dr J r 2 rm 
The minimum condition then implies that for r c 3> £ 
1 dfl, 







2ix da 



a=l 



m 2 p m (r c y - I dr-(l - pi) 2 ~ m 2 - E 3 (m) , (39) 



i.e. 



E 3 (m) = m 2 . (40) 

Once (^3), ( ^8|) is numerically solved for a given m, the numbers Ei(m) and E 2 (m) can 
also be evaluated. Ref. |39j supplies a value £2(1) — 0.279. This value is confirmed by 
ref. O which furthermore supplies E 2 (2)/4 + E 3 {2)/8 = 0.604, i.e. E 2 {2) = 0.416 and 
£"1(1) = — In 2.24 = —0.806. The mentioned values of EAm) are summarized in table II. 
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3.3 Interaction of vortices at large distances 

Let us now calculate the interaction energy of a given vortex configuration defined by its 
potential $ fl2"4"D , taking the core structure of the vortices into account. The interesting 
term in the hamiltonian ( pT|) is the third one which couples p and 9. Inserting (p2[), we 
obtain 

Jd 2 rp 2 \V6\ 2 = Jd 2 rp 2 \V$-n x V$| 2 

= Jd 2 rp 2 {\V$\ 2 -2W- (n x V$) + |V$| 2 } (41) 
= J d 2 r {p 2 | W| 2 - 2tfn • (Vp 2 x V$) - $(p 2 V 2 $ + Vp 2 ■ V$)} . 

The boundary terms appearing in the partial integration which leads to the third line of 
( |4T| ) would yield a positive contribution to the total energy which diverges logarithmically 
with the system size unless the vortex configuration is neutral (J2i m i = 0), in which case 
they vanish. Non-neutral configurations are therefore thermodynamically suppressed (at 
least at low temperatures), and we consider here only neutral ones, for which ( |4~1~D is 
correct. 

To find the saddle point field configuration for given vortices we now have to minimize 
( pl| ) (with (pLTD inserted) with respect to p and Because of (|23|) and p(Ri) = the 
p 2 V 2 $ term vanishes identically. In general the other terms interact in a complicated 
manner. However, we are interested at vortices at distances > £ since we will treat the 
short-range fluctuations by other means in section 4. Let us for simplicity assume that 
the smallest distance |Rj — Rj| between any two vortices in the given configuration is 
R ^> £ and look for the leading terms in an expansion in C./R. The vortex cores are now 
well separated and in any of them p is expected to be given by the isotropic 1-vortex 
solution of subsection 3.2 plus a correction of the order 0{^/ R) (a justification will be 
given below). Together with (|24]) and by symmetry arguments this implies that Vp 2 x V$ 
is of order 0(^/R) which in turn gives (see (|4TD) a saddle point configuration of $ of order 
0(£/R). (0) therefore yields 

J d 2 rp 2 \V6\ 2 = -Jd 2 r $(V$ • Vp 2 ) + 0(£ 2 /R 2 ) 

= -e j ^ (r r R R.p p2 in + °^i r2 ^ ■ ( 42 ) 
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In the second line, simply ( p4|) has been inserted. The important contributions to the 
integral in fl42|) stem from a region of size O(^) around Rj, the rest is again of order 
0((, 2 /R 2 ). To calculate this integral, we may then assume that Rj = and replace p by 
the 1-vortex solution with vorticity rrij centered at the origin, p m . (see subsection 3.2). 
We further have to distinguish two cases: 



1. i ^ j; i.e., Rj = Rj — Rj =: AR, |AR| > R. In this case the logarithm in (|42|) 

is dominated by a constant term ln(|AR|/£) and contains further terms of order 

0(^1 R) which couple p in the core of vortex j to the other vortices. This leads 

to an a posteriori justification of the above assumption that the corrections to the 

1-vortex solution p m . around Rj are of the order 0(£/R). Again invoking symmetry 

arguments we finally obtain 

, 2 r-Vp 2 |r-AR| |AR| r°° dp 2 

err In- = 2ti In - I dr—^ + O^/R 2 ) 



|r| 2 £ £ Jo dr 

7 



I AT? I 

27rln^=P + C(e 2 /i? 2 ) , (43) 



since p m , (r) — > 1 for r — > oo. 



2. i = j; i.e., Rj = Rj = 0. This gives the strong "self-interaction" of the core 
amplitude of a vortex with its own phase configuration, 

2 r-Vp 2 |r| _ ^ r.Jp^.T , //Vt2/D2 , 



(i 2 r , ; ln^ = 2tt / In - + 0(C/R 

\r\ 2 £ Jo dr £ 



-27rS!(m) + 0(£ 2 /^ 2 ) • (44) 



We can now insert these results into fl4*2|), add the terms of ( pT|) which contain p only 
and express everything in terms of 1-vortex quantities. So if ipo = ^({Ri? m i}) is the 
minimum-energy field configuration with given vortices of vorticities rrii at positions Rj, 
its energy is given by 

#[Vo] = "4vr £ m 4 m 3 In < R ' ~ R ^ - 4vr £ p GL (m,) + 0( ? 2 /i? 2 ) , (45) 

i<j 

where 



2 

777- 

-2fi GL (m) = m 2 E 1 (m) + E 2 (m) + — . (46) 

To Oth order in £/R, this is just the hamiltonian of a neutral two- dimensional Coulomb 
gas of integer charges whose core energy (or chemical potential p GL ) depend on mj. 
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Note that in particular (45) contains two-body interactions only. We will later absorb the 
overall factor 4tt in the effective Coulomb gas temperature. 

For later use we also calculate the quantity £~ 2 / <i 2 r(l — \ipo\ 2 ) by relating it to H[ip Q } 
in the following way: since ipo is the minimum-energy field configuration for a given 
vortex configuration {Rj, mj} as boundary conditions, the functional derivative 5H / 5iPq(y) 
vanishes nowhere in space except at the vortex centers r = Rj. On the other hand, ipo{ r ) 
itself has zeros precisely at the vortex centers, such that (with ip = pe ie ): 



° = / rf2r ^ = / d2r {"^ p2(1 " p2) + |Vp|2 + p2|W|2 } 



(47) 



Now after adding £ 2 J d 2 r(l — \ipo\ 2 ) on both sides, the r.h.s. becomes almost identical 
to H[i[)q] apart from an additional factor of 2 in front of the E 3 integral (see Q3"7D). In 



analogy to ( f45|) we therefore obtain 

/^(l-l^ol 2 ) = -4vr g m imj In |R * ~ Rjl - 4tt g (/i GL (m t ) - \m 2 ) + O^/R 2 ) 
= H[^] + nJ2m 2 + 0(e/R 2 ). (48) 

i 

Identity (f48|) will be of central importance in section 4 for the mapping of the full problem 
including fluctuations on a Coulomb gas. 

At low temperatures T < T v , the partition sum (|15|) will be dominated by configura- 
tions made out of bound vortex pairs of opposite vorticity. Furthermore, higher vorticities 
| m | > 2 have a much lower statistical weight than |m| = 1 for pair distances R £ be- 
cause of their higher core energy and stronger binding (see (f45l) ). We therefore assume 
that they play no essential role in the vortex unbinding transition and restrict ourselves 
in the following to configurations with 77ij = ±1. 



3.4 The functional derivative SH/Svpo 

The field configuration ip Q which we studied in subsection 3.3 is no real saddle point of H 
since it is subject to the vortex boundary conditions defined by {Rj, mj. However, clearly 
5H/5ip (ri) = for r ^ {Rj}, so we expect something like SH/Si/iofti) oc 8(r — Rj) for r 
in the vicinity of Rj. Since the functional derivative 5H/5tpo determines the linear energy 
of fluctuations away from ip (see section 4), we will study it here in some detail. We 
consider again a neutral vortex configuration with minimum vortex distance R ^> £. Let 
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us look at a neighborhood of the position Rj of the zth vortex and consider the particular 
fluctuation 5ipo generated by moving vortex i by an infinitesimal vector 5Rj. ( p3|) then 
implies that for |r — Rj| <C £ (assume for simplicity that m > 0), 



^o(r) ~c m (^^j (49) 
with complex representations f , R { for r, Rj as in (p3|). Consequently, 

*0b(r) ~ (r - iVf ~* ^ + 0{£/R) . (50) 

Now (f4"5|) implies that up to terms of order 0(£ 2 / R 2 ) the total energy of the vortex system 
changes by an amount 

5H[^o] = -Fi ■ 5IU = - l - (F*5Ri + FrfR*) , (51) 

where 

F< :=4*m'Z™J toIw (52) 

is the total Coulomb force on the zth vortex, which is of the oder 0(1/R). On the other 
hand, we can of course write 

which after insertion of (|50D leads by comparison with ( |5TD to 

, * . m— 1 

Note that the factor in front of the 5 function is of the order 0(£/R), which will permit 
us to neglect the corresponding terms in the perturbation expansion of section 4. 

3.5 Interpolation to the XY model 

For convenience, we worked from the beginning of section 3 in the continuum limit a/£ — > 
0. However, when we include fluctuations around the saddle point configuration ipo in 
section 4, we will have to impose a finite lattice cutoff a. We will then want to expand 
around a suitable saddle point for finite a/£. For the fluctuation corrections of section 4 
we will take the discreteness into account mainly by a suitable cutoff in momentum space. 
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The question addressed in this subsection is, what are the necessary modifications of fl45[) 
respectively (48) for finite a/£? 

The asymptotic behavior of the vortex potential at large distances R^> £, a will not be 
affected by the discreteness of the system. However, one expects an appreciable effect on 
vortex core energies since in the core regions the field ip always varies on length scales down 
to the order of the lattice constant a. The extremest example is the XY or "plane rotator" 
model, which corresponds to the £/a — > limit of (||), (f|) where p = = 1. XY vortices 
are known to interact at distances 3> a by a potential 4tt ln(R/ a) and to have a well-defined 
finite core energy 47r/i XY (which is approximately given by half the creation energy of a 
vortex pair centered at neighboring plaquettes of the lattice), whereas ( f45[) would predict 
a logarithmic divergency of /^ XY with £/a. More precisely, the energy calculated from (|6]) 
respectively (^) of a neutral pair at distance R 3> a, £ has the asymptotic forms 

{/ R \ 
Aw ( In 2/i XY ) for ^Ca 
/ R \ . (55) 

An I In — — 2fi GL j for £ > a 

where the respective chemical potentials are given by 

- 2/i XY = 7 + — w 1.617 (see ref. [21]), 
-2/i GL = E 1 (l) + E 2 (l) + ^ -0.027. (56) 

Here 7 = 0.5772 is the Euler constant and -E'i(l), -^2(1) are given in table II. The factor 
in front of the logarithm in (^) is not affected by the lattice regularization. 

We now look for a suitable interpolation formula for V^(R) which correctly reproduces 
both limits (p5|). Since we require an asymptotic behavior V^(R) ~ 47rlni? + const, for 
all £, it should be of the form 

Vfc(r) ~ 4vr {in I - 2p XY - F(£ 2 /a 2 )} , (57) 
where F is an interpolating function with the asymptotic behavior 

F(0)=0, F(X)~llnX + 2Gu GL -// XY ) for X>1. (58) 
A simple realization of conditions ( 58 ) is 

F(X) := ~ In (l + Xe 4(/iGL ~ MxY) ) ss - In (1 + 26.78X) . (59) 
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We tried to improve this simple guess for the interpolating function F(X) by fitting the 
first derivative F'(0) to the "true" value estimated by other means. However, this kind 
of improvement had no substantial effect on our final results, so we disregard it here. 
Finally, using expressions (|57|), (|59|) , the "saddle point" field energy (|5|) can be rewritten 
for vorticities m 8 = ±1 in the interpolated form (for R 3> £, a) 

H[ifj ] » -4vr m i m i ln ~ Rj * - 27riV(2/i XY + F{X)) + 0{a 2 /R 2 ) , (60) 

i<j a 

where N is the total number of vortices and X = £ 2 /a 2 . 

Starting from the lattice analog of (|4"T| ) one can work out an equivalent interpolation 
formula for the integral ((48|) , which however leads only to negligible corrections. We will 
therefore later simply use the expression (f48|) with the interpolated H[ip ] inserted. 



4 Thermal fluctuations around the "saddle point" 

As we already argued in the introduction we will identify the "saddle point" field con- 
figuration ipo with the long-wavelength components of the field ip with wave numbers 
|k| < 7i /a. Note that the latter still contain long- wavelength phase fluctuations which are 
absent in the "saddle point" field configuration. However, as we will argue in subsection 
4.2 they are not important for the transition and can be disregarded. We will now calcu- 
late the effective action S[ipo] (equations (T3), (fl4])) in a Gaussian approximation, i.e. we 
expand H[ipo + ^i] with respect to the fluctuations ipi up to second order and treat their 
coupling to the vortices perturbatively. 



4.1 Perturbational treatment of fluctuations 

The first step is to expand the GL functional (|]) around ip : 

H[^o + H = /^^(l-l^o + ^l^ + lV^o + V^I 2 ! 

+ /^{2(^ 1+ ^i*)^i| 2 + ^i| 4 } (61) 
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where the terms are ordered in ascending order in tpi. In a vortex-free configuration we 
would have = 1 and the linear terms in ipi would vanish. Now write again tp = poe l6 °, 
take all terms which are quadratic in ipi (and ifi*), replace po by 1 and collect them in 
some "free" fluctuation hamiltonian Hq. ( |6~ID can now be rewritten as 

H [i{> + Vi] = #[Vo] + #c#i] + Hj[ Po , 9 , fa] (62) 

where 

H [fa] = Jd 2 r ( 2 |^| 2 + (e^Vi) 2 + (e' e 'Vij 2 ) + |V^| 2 } (63) 

is the free fluctuation hamiltonian and all remaining terms are collected in the "pertur- 
bation" 

H, M M = /A^ + jgtf} (64) 

- /0(1 - A) + (e-*^) 2 + (^O 2 } + °W 

which describes the coupling of fluctuations to vortex cores and vanishes in a vortex- 
free configuration. Closer examination of (|6~3D reveals that, locally, the field ip has one 
massive and one massless component, the orientation of the local reference frame being 
determined by the phase 8 Q of the "saddle point" configuration. More precisely, the 
real and imaginary parts of the "gauge-transformed" field e~ ie °ipi correspond to massive 
amplitude fluctuations (mass 2/£ 2 ) and massless phase fluctuations, respectively. In 2 
dimensions the real-space propagator of a massless field is infrared divergent, which causes 
serious problems in a perturbation expansion. In the next section, we will therefore be 
concerned with the infrared regularization of the phase fluctuations. 

4.2 Infrared regularization of phase fluctuations 

In order to understand the physical role of long-wavelength phase fluctuations, we return 
to the representation (|21|), (f4"I|) of H, where phase fluctuations are described by a field 
The terms which couple to i? in (fill) only contain wavelengths <£, so long-wavelength 
components i9k with |k| l/£ decouple from the vortices and therefore play no role in 
the vortex unbinding transition (note that in the continuum XY limit £ a, decouples 
completely!). This leads us to assume that all components with |k| < n/a are irrelevant 
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and can be disregarded. We will then treat the components with it /a < |k| < n/a 
perturbatively. For technical reasons it is convenient to realize this infrared momentum 
cutoff ir/a in a "soft" way, by assigning a mass to the phase fluctuations in ([k|). For 
simplicity we choose this mass to be the same (= 2/£ 2 ) as for the amplitude fluctuations, 
such that (^) is replaced by the simpler, translationally invariant expression with isotropic 
mass matrix 

H [H= /d 2 r(-||^i| 2 + lWi| 2 l (65) 



v£ 21 

(We will argue in subsection 4.4 that the relation between "hard" and "soft" cutoffs is 
given by a 2 = a 2 + 27r£ 2 ). With this modification, the hamiltonian (|62"D now describes 
a system of interacting vortices together with a massive, harmonic field ipi (made up 
of amplitude and phase fluctuations and described by (|65"D ) which is scattered by the 
vortex cores. In the calculation of the effective action ([14]) we may now integrate over 
all components ipi^ with |k| < n/a since the infrared cutoff is taken into account by the 
mass. 

The free propagator corresponding to hamiltonian ( |65"D is given by 



where 

d 2 k 



Go(r) = / j^G (ky kr , (67) 



(27T) 



<W = • (68) 

For later use in the perturbational expansion, we derive here some properties of Go- The 
local mean square value of ip\ is closely related to the fluctuations already calculated in 

f < W |> = G (r = 0) = / ^ „ _L ln(1 + 2rf ) , (69) 
where we again used the notation X = /a 2 . More generally, one can show that 

— ln(l + 27rX) for r = 

Gu(v)^l -— InM for a < |r| < ^ , (70) 
2tt ^ 

for |r| > a, ^ 
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where the second line is relevant only if a < £. The n-fold convolution of Gq with equal 



arguments at the ends can be calculated in analogy to ( 69|) , 



G» +1 ) (r = 0) = / ■ ■ • / G (l, 2)G (2, 3) ■ • ■ G (n + 1, 1) 

f ^g (kv^ ~ -i r Va2 d(fc2) 

7 (2tt) 2 01 ; ^47r7o (2/C 2 + fc 2 )™ +1 
1 f/2\ -n /4tt 2\~ 



(71) 



Ann \\£ 2 J \a 2 £ 2 
which is never larger than its £ 3> a limit: 

f (r = 0) < . (72) 



e n 



We write out the particular case n = 1 explicitly: 



rf 2 rG (r) 2 « ^ — < £ . (73) 

° V ; 4 1 + 2ttX - 8tt K ' 



4.3 Diagrammatics 



We will now calculate the effective action S[ip } (eq. (0)) by diagrammatic perturbation 
theory, starting from the decomposition (|B2"D, Q61D, flB5| ) of the original hamiltonian. The 
propagator of the theory is yGo, defined by (|67D, (|68|) . Furthermore, up to order OijpD 
the interaction term — yiij is represented by a sum of the diagram elements shown in 
fig. Bl with combinatoric factors of 1/2 for the symmetric diagrams (d) and (e). The 
corresponding analytic expressions read 

(a) = -Kfgi; 
K J 2 J <tyo 

^ K f 5H I 

(b) = -jjw^ 1 

(c) = ^[(i-pI)\H 2 



e 

(d) = ^ / (1 - Pl)e 2ie ^f (7.1) 



2£ 

(e) = £ / (1 - f%)e-**>tf 



2£ 2 

respectively. Omission of terms of higher order than Ofyf) just means that we are working 
in the gaussian approximation. The perturbation series for the effective action S[ipo] (eq. 



(|T4])) may now be expressed as — yif^o] plus the sum of the diagrams shown in fig. 
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with the appropriate combinatorial factors (plus higher order terms to be discussed in a 
moment). We calculate these different diagrams separately, keeping only the "leading" 
terms, i.e. those which contribute to the logarithmic interaction or the chemical potential 
in H[i[)q\ (eq. (|45|)), and omitting terms of order 0(£ 2 /R 2 )- 
Consider first the "chain" diagrams in fig. [|: 

« jG (r = 0)(J^\\F i \* = O{?/#), (75) 

where we used ([54]) for m = 1. Similar arguments show that all of the "chain" diagrams 
((b), (c), (d) and "longer" ones) are likewise of order (9(£ 2 / ' R 2 ) and they will therefore be 
neglected. 

The "loop" diagrams (e), (f), ... in fig. ^ are more interesting. The only one which 
gives a contribution to the logarithmic term in the vortex interaction is diagram (e): 



(e) = 2G (r = 0)r 2 /(l 



Pi) 



ln(l + 2nX)(H[^ ] + nN) + 0(a 2 /R 2 ) (76) 

Z7T 



where we used fl69|), fl4q) for rrii = ±1 and N is the total number of vortices in ipQ. The 
higher order loops all give positive contributions to the vortex chemical potential. The 
most important of them is 

^x(f) = 1 I /g (1,2) 2 (1- Po (1) 2 )(1-Po(2) 2 ) 



2 e 4 AJ 2 

< 2 r 2 (/ gi) r 2 /(i - pD 2 - ^i^f + ^ 2 /R 2 ) , (77) 

where the < sign expresses that we used a Schwarz inequality to estimate the integral 
at the l.h.s., and we used (|73|), (|37|), (|40| ) to evaluate the r.h.s. (note that we took 
the symmetry factor of the diagram already into account). Similarly we treat the next 
diagram, 

\x(g) = ^ I /^Go(l,2) 2 (l-p (l) 2 )(l-p (2) 2 )cos2(^ (l)-^ (2)) 

^ -TT^x +0(a /R ) < (78) 

where in addition we neglected the cosine factor on the l.h.s., which is of course a more 
serious approximation; however, the contribution of ( 78 ) to is smaller than that of ( 77Q 
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by a factor of 1/8 in any case, so we neglect this error. We note that both approximations 
(|77|), (ff8|) would be good in a limit where the range of Go( r ) 2 is much smaller than that 
of (1 — po( r ) 2 ), but (|70"|), fl32| ) imply that both ranges are actually of the order of £. 

The dominant class of terms in the whole series is given by loop diagrams of the form 
shown in fig. [| with n + 1 nodes, n = 1, 2, and a symmetry factor of l/(n+ 1) (the first 
two are (f) and (h) in fig. |]). An upper limit to the contribution of any of these diagrams 
may be found by again using a Schwarz inequality as above, and directly applying (|72|): 

1 x (Fig. 5) = 



n + 1 

1 / 9 \ n ^ 

+ 1 ( e ) / " l+i Go(1 ' 2) ' ' ' Go( " + !. !)<! - ^W 2 ) • • • I 1 - /*>(" + !) a 



< 



27in(n + 

The integrals on the r.h.s. monotonically tend to zero with increasing n and their coeffi- 
cients alone form a rapidly converging series with sum 1/27T (for comparison, the first term 
corresponding to diagram (f) in fig. f| already contributes l/47r). Since we already overes- 
timated the contributions of (f), (g) we therefore neglect all higher "loop" diagrams and 
hope for the best. We will see in the next subsection that any further positive contribution 
to the vortex chemical potential amplifies the observed effect and thereby strengthens our 
argument anyway. Putting everything together, we now obtain 

*o] = -{j-^ln(l + 2 J rA')}ffW + jln(l + 2 J rA') 

which together with ( |60"1) gives the effective action of the interacting vortex system. 

We finally note that of course all loop diagrams in figs. |], |5| can be calculated numer- 
ically to an arbitrary precision, but in view of the many approximations and qualitative 
arguments involved in our derivation of S we considered the above rough estimation to 
be more adequate. 

4.4 Vortex phase space division and effective Coulomb gas pa- 
rameters 



The result (^) of the last subsection already correctly defines the partition sum (|T5| ) of 
the vortex system. However, the unit of length in ( [To] ) is a (which corresponds to R c in 
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eq. fllTD)) so that in order to make the link to treatments of the neutral Coulomb gas in 
the literature we have to write the effective action in the form 



CG 



^ | R . pi 

S[M = H m i m 3 111 % ~ ^ + N 111 Z 

whereas in (|60|) the original lattice spacing a appears in the logarithm. (|8lD then defines 
the effective dimensionless Coulomb gas temperature T CG and fugacity z CG . Using the 
neutrality condition J2i m i — and the fact that m; = ±1 we can rewrite the sum in (|60|) 

as 

\ - , |Rj — Rj| x - , |Rj — Rj| N a 

} j mitrij In — = 2^ mirrij In — - In — . (82) 

The second term will contribute to the fugacity z CG , so we have to express a as a function 
of the model parameters a, £. Recall that a was an infrared cutoff on phase fluctuations 
defined so as to have the same effect as a mass 2/£ 2 . To determine a from this condition, 
note that the most important term in the perturbation series is (f) in fig. [| (and eq. ( [74]) ) 
since it is the only one which contributes to the logarithmic vortex potential. Its size is 
mainly determined by the factor 

f<w 2 >=iH i+2 ^) (s3) 

calculated in (^). Now if one had chosen the "hard" momentum cutoff ir/a instead of 
the mass, the local fluctuations would have been given by 



K 2 r d 2 k 1 1 r^/" 2 d(k 2 ) 1 



7,2 



l\<k x , y <l (2vr) 2 k 2 W 4vr .L/a 2 k 2 4vr ln a 2 ' (M) 
Identifying this expression with (|83| ) leads to 

a 2 \ a 

or, equivalently, 



ln- = -ln(l + 2vr^ (85) 



a 2 = a 2 + 2vr£ 2 . (86) 

We insert (|^) successively in (0), (0), (|80| ) and obtain the required form (|8lD of the 
effective action, the Coulomb gas parameters now being given by the following expressions: 

I 



T 



CG 



27r J FT-21n(l + 27rX) , (87) 



ln^ CG = ^{2 /UxY + F(X)-iln(l + 2vrX)} 



1 . . 9tt X 

+-ln(l + 2vrX) + — - 



8 1 + 2ttX 
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where 2/z XY = —1.617 (see eq. (|56D ), X = ^ 2 /a 2 , and F(X) is denned by eq. fl59|) . 

Equations fl87|), fl88|) are the final results of this paper. The expression in curly brackets 
may be interpreted as a renormalized vortex chemical potential, whereas the second line 
of (jS8|) contributes to the renormalized phase space division. 

In fig. |6| the z CG vs. T CG curves are drawn in Minnhagen's Coulomb gas phase diagram 
for different values of the GL correlation length £. With increasing £/a the curves cross 
the phase boundary at increasingly higher values of z CG , finally reaching the first order 
regime for values of £>gl Since we argued in subsection 2.1 that in superconducting 
films always £ > a, we conclude that the latter are good candidates for a first order vortex 
unbinding transition. Because of the approximations and relatively rough estimations 
involved in our calculations we do not insist in the quantitative information conveyed by 
fig. |6|. However, we believe that the trend is clear enough to be reliable. 

Another interesting (although not unexpected) result might be that for increasing 
£/a the transition occurs at increasingly higher values of K (which according to (Q) 
corresponds to smaller values of the "physical" temperature) . In other words, the distance 
between T v and T c o increases with increasing £ /a. 

5 Summary and conclusions 

In this paper, we investigated the nature of the transition in the two- dimensional Ginzburg- 
Landau model of a neutral superfluid. In doing this, we assumed — as is usually done 
- that the transition is caused by the collective unbinding of vortex pairs. However, 
we paid particular attention to short-wavelength amplitude and phase fluctuations of the 
order parameter on length scales <£ (the Ginzburg-Landau correlation length), which 
we showed to be strongly coupled to the vortices. This is in contrast with previous work 
in the literature where amplitude fluctuations usually are neglected altogether and phase 
fluctuations are assumed to be only weakly coupled to the vortices. 

Eliminating perturbatively these short- wavelength fluctuations, we derived an effective 
free energy for the vortex degrees of freedom. We argued that this effective vortex gas 
still is a 2D Coulomb gas (i.e., the vortex- vortex interaction varies logarithmically at large 
distances), however both the effective temperature and the fugacity of the Coulomb gas 
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are strongly renormalized if £ is larger than a microscopic cutoff length a for fluctuations. 
We argued that these considerations should be relevant for those superconducting films 
for which a BCS description is qualitatively correct, since BCS theory predicts that always 
£>a. 

By this elimination process of small-scale fluctuations we furthermore obtained a mi- 
croscopic interpretation of the effective vortex phase space division A as an interesting 
secondary result. In particular, we could clarify the relation between A and the Ginzburg- 
Landau correlation length £. 

Our subsequent conclusions concerning a possible first order transition then relied com- 
pletely on the correctness of Minnhagen's self-consistent theory of a dense 2D Coulomb 
gas, which predicts a first order vortex unbinding transition with non-universal (not KT 
like) properties for sufficiently large values of the vortex fugacity (z CG > 0.05). Unfortu- 
nately, the rather subtle differences between a KT like and a first order vortex unbinding 
transition seem to be outside the scope of present experimental investigations of real su- 
perconducting films. However, we think that an interesting and probably feasible (though 
likewise very hard) problem would be to numerically investigate the critical properties of 
a lattice Ginzburg-Landau model including fluctuating order parameter amplitudes. 
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Tables 



model 


R c 


A 


/i CG 


XY 


a 


a 2 


-0.809 


GLCG 


2.24£ 


16.4£ 2 


-0.390 



Table I: Effective CG parameters for the XY model (values from ref. pif ) and for the 
"bare" GLCG (values from refs. Q §5J). 



m 


Ex 


E 2 


E 3 


1 


-0.806 


0.279 


1 


2 




0.416 


4 



Table II: Some numerical values for the E parameters of the 1-vortex problem (eqs. fl35] 
(|36|), (|37|)), taken from refs. |39|, pD| . See explanation in the text. 
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Figure Captions 



Figure 1: Plot of Minnhagen's generic CG phase diagram 32, 33]. Also shown are 
the z CG (T CG ) relations (eq. fli"8D ) for the XY and GLCG models, using the parameters of 
table I. In both KT like vortex unbinding transition is predicted. 



Figure 2: Same as fig. [I], but with several different values of the cutoff R c of the vortex- 
vortex interaction (see eq. fllTD). The qualitative properties of the vortex unbinding 
transition of the XY and GLCG models are not changed. 



Figure 3: Diagrams corresponding to the interaction term Hj of the hamiltonian, eq. 
Figure 4: Diagrams corresponding to the gaussian approximation of the effective action 

Figure 5: Class of diagrams represented in eq. (|79|) . 

Figure 6: Again Minnhagen's CG phase diagram as in figs. p], but now with our final 
z CG (T CG ) relation for the GLCG (eqs. ©), plotted for different values of (/a. The 
limit £/a = corresponds to the XY model as shown in fig. |]. With increasing £/a the 
z CG (T CG ) line is shifted towards the upper part of the phase diagram, finally reaching the 
first order part of the vortex unbinding transition. 
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